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Preface 


PRECEDING PAGE EIA:'^K NOT FILMED 


This technical memorandum has been prompted by recent changes to the way 
in which solar radiation is modelled in the g'=*neral circulation model (GCM) at 
the Goddard Space Flight Center's Laboratory for Atmospheric Sciences (GLAS)* 
The revised solar radiation treatment combines a direct-diffuse radiative 
transfer model, developed by the author, with much of the orginal model, 
developed by A. A. Lads and J» £• Hansen. It is a pleasure to acknowledge the 
substantial contribution to the present form of the code that has been made by 

Dr 8. Lads and Hansen. Much of their original code remains unaltered. The 

author, however, accepts sole responsibility for any errors or omissions that 
may be present in this documentation. 

I would like to express my appreciation to Dr. M. Haiem and Dr. J. Shukla 

for providing the encouragement and opportunity to complete this work. 


R.D. 
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1. INTRODUCTION 


The differential absorption of solar radiation by the earth-atiaosphere 
system initiates various dynamical and physical processes responsible for 
maintaining the general circulation of the atmosphere as well as local climatic 
regimes. In any numerical model which attempts to directly simulate these 
processes, it is therefore necessary to include as a sunset of the model a 
parameterization scheme for the absorption of solar radiation which is both 
accurate and efficient. 

The general circulation model (GCM) in use at the Goddard Laboratory for 
Atmospheric Sciences (GLAS) generates most of the interactive atmospheric 
variables that affect solar radiation; namely cloud cover, humidity, temper- 
ature and pressure. These variables are supplemented in the parameterization 
scheme described here with climatological values of ozone and surface albedo 
to yield absorption values for the surface and nine atmospheric layers on a 
4° latitude by 5° longitude horizontal grid. 

The absorbed solar radiation is then incorporated into the d^abatic terms 
of the next time step of the GCM, and a new round of atmospheric variables 
are generated. In practice, the solar radiation routines are called every 
third time step, corresponding to 30 mins, of model time, in the Interests of 
computational efficiency. 

This memorandum documents the parameterization as developed by mid 1981 . 
The development rests heavily on the work of Lads and Hansen (1974), with 
extensions to account for solar zenith angle effects on the earth surface 
and a revised radiative transfer treatment in clouds. 

Briefly, the current parameterization allows for ozone, water vapor and 
surface absorption under a variety of conditions and for arbitrary sun angle. 
Rayleigh scattering by the clear atmosphere is parameterized and multiple 



scattering by cloud droplets calculated layer by layer. The spectral Inte- 
gration for regions with water vapor absorption is reduced to a summation over 
five terms using a probability distribution. Limitations to the current 
scheme include ignoring the effect of aerosols, absorption by cloud droplets 
and the spectral nature of the surface albedo. 

In the description that follows, the scientific background for the 
radiation calculations is first summarized, followed by a description of the 
FORTRAN code. Comparative results between the old and new parameterizations 
are then presented. 


2. SCIENTIFIC OVERVIEW 

This section discusses the scientific basis on which «:he radiation routines 
have been developed. Much of the original treatment of Lacis and Hansen (1974) 
remains unaltered, and is therefore simply summarized. Developments to the 
parameterization scheme since the Lacis and Hansen version are discussed fully. 


2.1 Input Variables 

Before the transfer of solar radiation can be modelled, a number of input 
variables describing the extraterrestrial insolation and the radiatively active 
components of the earth-atmosphere system must first be defined. The usual 
input variables in the GLAS GCM, with their FORTRAN counterparts in parentheses, 
are: 


(LAT) 
d, JDAY 
th, (TOFDAY) 

(AG) 

(CL) 


- latitude in degrees. 

- time of year, in Julian days. 

- local apparent time, either specified explicitly or 
obtained from GMT, longitude and the Equation of Time. 

- conventional surface albedo 

« type of cloud, if any, present at each level of the atmosphere. 





q, (SHL) - specific humidity at each level of the atmosphere. 


T, (TL) 


- mean temperature (K) for each level of the atmosphere. 


Pg, (PRESS) - surface pressure, (mb). 

The following additional variables may then be readily derived: 

Mo> (COSZ) - cosine of the solar zenith angle. See section 2.2. 

UqSo* (SCOSZ)- the extraterrestial irradiance (ly/day) = See section 2.2. 

M, (COSMAG) - a magnification factor to account for the slant path and 
refraction of the direct beam. Given by 
M - 35/(1224 mo^ + 

(OZALE) - the ozone amount in cm (NTP) above each atmospheric level. 
Determined from climatological values using subroutine 
OZONE2, with JDAY and LAX as Independent variables. 

Pjj, (PLE) - pressure (mb) at the top of the nth atmospheric level, 
determined by 

Pn“Pl + (n-l)(Pg-Pj^)/N, where, typically, Pi«10mb and N*9. 

yjj, (SWALE) - the effective water vapor amount above each atmospheric 
level. See section 2.4. 

Ag, (RSURF) - the solar zenith angle dependent component of the surface 
albedo. See section 2.5. 

Tq, (TAUC) - the cloud optical thickness in each atmospheric level . 

See section 2.6. 


2.2 Astronomy 

The extraterrestrial solar irradiance, Sq, and the solar declination angle, 
6, vary slowly with the time of the year. Following Paltridge and Platt (1976), 


Ch. 3, for example, the time of the year may be expressed in radians as 

'*'n * 2ird„/365 2.2.1 

where d^ ranges from 0 on Jan 1 to 364 on Dec 31. 

An approximate expression for 5, accurate to 0.0006 radians is then 
6 - 0.006918 - 0.399912 cosY + 0.070257 sinY 


- 0.006758 cosZY + 0.000907 sin24' 


2 . 2.2 


- 0.002697 cos3i + 0.001480 sin3'i' 


•■■il 
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The extraterrestrial solar Ir radiance when the earth Is at a distance R* 
from the sun is given by 

So « 2823.9 (T[*/R*)^t 2.2.3 

where R^ is the mean earth-su« distance and the solar constant has been taken 
ati ;/tp fr >m Willson et al. (1981). The factor (R*/R*)^ in eq. 2.2.3 may 

be deU rr.luad to an accuracy better than 10“^ by 

(R*/R*)2 - 1.000110 + 0.034221 cosf + 0.001280 sin^ 

+ 0.000719 cos 2Y + 0.000077 sln2‘i' 2.2.4 


A further astronomical variable of interest is the Equation of Time, Eq, 
given in radians by 


Eq = 0.000075 + 0.001868 cos'P « 0.032077 sin^ 
-- 0.014615 cos2^ - 0.040849 sin24^. 


2.2.5 


from which local apparent time is given by 

local apparent time = local mean time + Eq. 2.2.6 

Given the local apparent time t^, in radians, measured from local apparent 
noon, the solar zenith angle 0 q for latitude <p is given by 

cosOq = sin6 sin(p + cos6 cos^ cos tj^. 2.2.7 


2.3 Absorption by Ozone 

The absorption of solar radiation by ozone closely follows the original 
treatment of lacis and Hansen (1974), making use of the assumptions that 

(a) ozone absorption takes place above any levels of significant scattering. 

(b) ozone absorption and water vapor absorption occur in separable spectral 
regions. 

Absorption is then a function of the effective path x, traversed by the 
radiation through the ozone layers. If x is in cm (NTP) , this function is given 
empirically by 

(X) =■ 0.02118/(1 + 0.042X+ 0-000323x^) 

O 2* 

+ 1.082/(1 + 133.6x)^-®‘^^ 


+ 0.0658/[l + (103.bX)^J 


2.3.1 


where the first term on the r.h.s. refers to absorption of visible radiation 
by the Chappuis band, and the other terms to absorption of ultraviolet radiation 
by the Hartley and Huggins bands. 

The effective path traversed by the direct beam to reach level n is X^BMUn 
(section 2.1), and the effective path traversed by radiation reflected from 
the atmosphere beneath the ozone layers is +1.9 (U^ - Un), where the 

factor of 1.9 accounts for diffuse transmission, and Ut refers to the total 
ozone amount. The total absorbed radiation in each layer containing ozone is 
then given from eq. 2.3.1 by 

Aq “ - AQ2(^n) HtAQ2(^n*^ ~ Aqz(^+1 2.3.2 

where R refers to the albedo of the underlying atmosphere to visible radiation 
and is determined for a cloudy atmosphere as '<n section 2.8. For cloud free 
conditions, however, an empirical formula for R, chosen explicitly to be 
consistent with eq. 2.3.1 is used as follows: 

R “ Rl + (1 -Ri)( 1-R2) Ag/(l-AgR2) 2.3.3 

where Rj = 0. 2186/(l+0.816no) 
and R2 = 0.144. 


2.4 Absorption by Water Vapor 

The fractional absorption for solar radiation traversing an effective 
water vapor amount x (centimeters of precipitable water vapor) was found 
empirically by Lads and Hansen (1974) to be given by 

K.v.M ~ 2.9/[(l+141.5x)°*^^^ + 5.925x] 2.4.1 


The effective water vapor amount above each level is found by Integrating the 
specific humidity with respect to pressure, making a rough scaling correction 


for the effects of temperature and pressure on absorption. That is, 
yn - - ^ q(P/1013)(273/T)^ dP 

go 


2.4.2 



2.4.3 


I 


The direct solar beam then traverses an effective water vapor amount 

to reach the top of layer n, and radiation reflected from the surface traverses 
an amount 

’tn+l* • MXj, + 5/3 (y^-yn+p 
to reach the bottom of layer n. 

If we Ignore scattering* the total absorption due to water vapor in each 
layer is given by 

^n “ I^oSqC A v.v. (*^+l) “ •^gt^.v.^^n ^ ” Ay,.v.(Xn+l )3^ 

2.4.5 

Eq. 2.4.5 is used only when the atmosphere is cloud-free. For cloudy atmosplieres 

there is usually significant scattering of solar radiation at those wavelengths 

where water vapor absorption is also significant, and a different approach is 

used, with eq. 2.4.1 being replaced by 

A„.v.(x) = 1 - p(ki)e“V 2.4.6 

i=l 

where p(k)dk is the fraction of incident flux associated with an absorption 
coefficient between k and k+dk. The values of p(kj^) and in the current 
version of the GI-AS GCM are given in Table 2.4.1. 

Table 2.4.1 

Discrete probability distribution of water vapor 
absorption coefficients 


i I ki 
I 

1 


P(ki) 


1 

1 0.005 

' 0.107 

2 

' 0.041 

, 0.104 

3 

1 0.416 

' 0.073 

4 

1 4 32 

1 0.044 

5 

! 72.' -1 
1 

‘ 0.025 

1 i 




The absorbing optical thickness due to water vapor In layer n Is then given 
for each 1-1,2,. ..,5 by T„.v.,n,l * U(yn+1 " 2.4.7 

The absorption optical thickness due to water vapor, Ty.v.» scattering 

optical thickness due to clouds, are then combined to give the total 
optical thickness for each layer and each value of 1. That is. 

"^n,! “ '^w.v.,n,l + ^c,n, 2.4.8 

with the effective single scatter albedo for the layer being given by 

an,l • Tc,n/Tn,l 2.4.9 

and the subsequent radiative transfer problem Is solved as discussed in sections 
2 . 7 and 2.8. 

2.5 Absorption by the Earth's Surface 

The surface albedo Is divided Into two components, one of which, Ag, Is 
applied to the diffusely Incident radiation, and the other, Ag(pp) , is applied 
only to the direct solar beam. The two components are re ated (see for 
Instance Paltrldge wd Platt 1976) by 

Ag(0o) * Ag + (1-Ag) exp [-0. l(9O°-0o) ] 2.5.1 

For cloudy skies, the surface Is treated as If it were a non-transmitting 
atmospheric layer In sections 2.7 and 2.8, and for cloud free conditions 
separation of the scattering and absorbing portions of the spectrum is again 
assumed so that 

surface absorption - PoSotO*647( l-Rj^y) + [0.353 - A^,.v.(x^j) ][l-Ag(uo) ]} 

2.5.2 

where the first term on the r.h.s. refers to absorption of radiation In the 
scattering (by Rayleigh atmosphere) portion of the spectrum and the second term 
to absorption of radiation in the same spectral regions as water vapor 
absorption. The combined albedo of the surface and Rayleigh a mosphere, Rq^v* 

Is given empirically by 

R, . 

(l-RjAg) 


2.5.3 



where 

Rl - 0.433/(1 + 6.43}Jo) 

2.5.4 

and 

R2 - 0.093. 

2.5.5 


2.6 Intensive Cloud Properties 

The GLAS GCM currently genisrates four types of clouds, namely: 
"supersaturation'* clouds in any layer except layer 1, "low level convective" 
clouds in layer 7 o- B, "mid level convective" clouds in layer 5 or layer 6, 
and "penetrating convective" clouds in any rour consecutive layers except 
layer 1 or layer 2. The only difference in the intensive radiative properties 
of these cloud types Is in their optical thickness, Tq, as given in Table 2.6*1* 


Table 2.6.1 



Optical 

Thickness of Different 

Cloud Types 

r 

1 

Cloud Type 

1 

1 C^oud Layer 

j _ _ 

1 

1 

J 


f 

1 

Supersaturation 

1 2, 3, ,.5. 6, 7, 3 or 9 

I 

1 

1 

1,2, 4, 6, 6, 8, 8, 8 | 

1 


1 

1 

(respectively) j 

1 

Low Level Convective 

' 7 or 8 

1 

1 

4 

16 

i 

1 

t 

Mid level Convective 

1 

1 5 or 6 

1 

4 

8 

1 

1 

Penetrating 

* 4 consecutive layers 

i 

1 

8 in each of the j 

1 

Convective 

1 between 3 and 9 

1 

four layers , 

u 


! 

u 

1 


The other int/'nsive properties, required for the radiative transfer calculations 
of section 2.7, namely the single scatter albedo and the asymmetry factor 
of the phase function for scattering by the cloud droplets, are assumed to be 
constant for all cloud types. Cur rent values in the GLAS GCM are = l and 

• B3 . 
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2.7 Radiative Properties of a Single Laver 

This section describes how the reflection and transmission characteristics 
of a single atmospheric layer are determined, given as Inputs the total optical 
thickness t, the single scatter albedo a, the asymmetry factor g, and the solar 
zenith angle convenience, we may also Introduce an extinction 

coefficient k and a vertical depth Zq to the layer, so that T»kZQ. 

The transfer of solar radiation through a single homogeneous layer is 
shown schematically In Fig. 2.7.1. The direct beam gives rise to diffuse 
Irradlance reflected upwards at z*>0, UO, and transmitted downwards at z=Zq, DO, 
as well as a depleted direct beam of Irradlance 

Mo^o' “ UqSo exp [-T sec Qq] 2.7.1 

For X ^ 8, UO and DO are found using the Delta Eddington approximation of 
Joseph tt oJL (section 2*7*1)» with the Coakley and Chylek version of the two 
stream approximation being used for x < 8 (section 2.7,2). 

The transmittance, TL, and the reflectance, RL, of the layer to diffusely 
incident radiation is also required, and here the Sagan and Pollack (1967) 
version of the two stream approximation is retained. 


Thus for a < 1 

RL - (u2-l)(l-e~^^'^)/den 2.7.2 

and TL - Aue'^'^/den 2.7.3 

where X » {3(l-a)(l-0g)}**, 2. 7. A 

U - (l-OS)‘*(l-a)“**, 2.7.5 

and den - (u + 1)2 - (u - l)^e“^^^ 2.7.6 

For the conservative case of a * 1, we have 

TL - U + |\l-g)t}"^ 2.7.7 

and 

RL « 1 - TL 2.7.8 
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fig. 2.7.1 
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2.7.1 Delta-Eddington approxlMtlon 

For T ^8 the Delta-Eddington approxiBation of Joseph tt aJt (1976) 
appears superior to alternative siaplf approximations and has been adopted 
for use in the GCM. 

The intensive parameters are first scaled so that 


k’ - k(l - af) 

a‘ - a(l - f)/(l - At) 

8’ • (g - f)/(l - f) 
where f«g^. 

The resulting solutions for UO and DO are then given by 
UO - 2it Io(0) So, 

and 

DO - 2n Io(2o> Sq. 

wnere Io(0) - cxp(>nZo) + + o^, 

lo(Zo) - Oj + exp(-nZo) + exp(-k*Zo/p(j) 

and 

n ■ {Jk'2(i - a’g'Hl - a')}^ 
uj • {UjM »2 - U|Wj exp(-nZo)}/{det) 

02 - - U1W2 exp(-nZo)}/{det) 

03 - -3k’2 a'[l + g'(l - a‘) j/4n(k'2/y^2 _ ^2) 

and 

det - exp(-2nZo) 

ui • - n 

vi ■ h + n 

wi “ “(/i + k’/uo)<»3 ~ 3k'wo“*9 

u >2 ■ -[(/i -• k'/po )“3 ■ 3k*Uo^'g'/4v] exp(-k ’Z q/mo) 

and 

h - 3k'(l - tt'g')/2 . 


2.7.9 

2.7.10 

2.7.11 


2.7.12 


2.7.13 

2.7.14 

2.7.15 


2.7.1b 


2.7.17 


2.7.18 


2.7.19 


2.7.20 

2.7.21 

2.7.22 


2.7.23 


2.7.24 

2.7.25 


-n- 
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2.7.2 The Two Stream Approximation 

For small optical thicknesses, the Delta-Eddington approximation is less 
precise than the two stream approximation of Coakley and Chylek (1975). which 
tends to the single scatter approximation as t-K). 

A generalllsed form of the two stream approximation leads to 



l» - So e"^^ + 02 ’ + 03 "} 

2.7.26 

and 

DO • Sq + <* 2 ^ e e 

2.7.27 

vliere 

1 «p(l - a)**(l - a + 206 i)** 

2.7.28 


0 ^- - n-/(pQ - X ) 

2.7.29 


■ {o 3 (Ci“C 2 )e ^^-030616 ^^*^°}/{Cj+C,-(Cj-C 2 )e } 

2.7.30 



2.7.31 


+ + + -Xt 

“2 *“"3 -“l ® 

2.7.32 


“2 “ “2 ^‘^l ” 

2.7.33 


O 0 J + (1 - 62 ^ 5 ? ■*'Po 

2 . 7 . 34 


n~ * -ai>o 8 j + 62 ^ (1 ■ ~Po 

2.7.35 


Cj * 1 - a + oSj 

2.7.36 

and 

C 2 ” X/p 

2.7.37 


iiie Coauiey-Chylek approacn nas been further modified to allow for 
an optimal choice of the backscattered fraction based on Wiscombe and 
Grams (1976), anu to make a delta-scaling of i, a and g to t', a' and 
g' analogous to the Delfa-Eddington scaling. 

Thus for X ^ 1 

P * 1/Po 2.7.18 

ei - P2(i^o) 
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and for 


T > 1 


P - 2 


2.7.40 


3 ^ - 1 / 2(1 


2.7.41 


wher ‘ Is chosen nuaerically from Fig. 3 of Wlscomhe and Grans (1976). 


2.8 M ulti-Layer Solution 

At each grid point of the GCM, the following quantities are first determined 

for each layer, n, taken In Isolation, using the methods of section 2.7. 

SO - the direct solar ir radiance through level n. 
n 

UO - the upwelllng Irradlance through level n due to SO . 
n n 

DO - the doimvelling diffuse irradlance through level n due to SO .. 
n n— 1 

- the reflectivity of layer n to unit diffuse irradlance. 

- ty transmissivity of layer n to unit diffuse irradlance. 

These quantities are then combined in the following scheme to give the total 

upwelllng and dovnvelling irradiances» U and D , through each level* from which 

n n 


the absorption in each layer is given by 


'n ■ “o -<Vl - Vl> 


n « 1, 2, ...» N 2.8.1 


wh^re 


.. is the total number of atmospheric layers. 

The scheme used is similar to that of Grant and Hunt (1968), with UO , 

n 

DO .. acting a' internal sources* and with some simplification due to the fact 

iH*i 

that the .ladiation is treated as though it were diffuse* 


We define D1 and U1 to be* respectl ;elv* the downwclling and up- 
n n 

welling diffuse Irradlances due to all radiation crossing level n which has 


J 





r- • 


oi- rc.:~ cuALrfY 


not previously crossed level n + 1. D1 and U1 may then be obtained 

n n 

recursively from 


D1 «(D1 , + UO ,CR J M ,T , + DO 

n r -1 n-1 n-1 n-1 n-1 n 


n = 2, . . . , N 


2 . 8.2 


U1 =(D1 R + UO ) M 
n n n an 


n = 1, . . . , N 


2.8.3 


where CR^ is the composite reflectivity of diffuse upwelling irradiance 
at level n due to ail layers above level n, and is given recursively by 


CR = R + T ,CR , T ,M , 
n n-1 n-1 n-1 n-1 n-1 


n = 2, .... N 


2.8.4 


and is the magnification factor for multiple reflections at level n 

between layer n and all layers above level n. N is given by 

n 

M /(I - CR R ) 
n n n 


2.8.5 


Note that 


CR^ = 0, DOj^ * 0, and Dl^ = 0. 


Once D1 and U1 have been found for all n, the total irradiances crossing 

each level, D and U , are found recursively be evaluating 
n n 


U = U1 + U M 

n n nrH n n 


n = N-1, N-2, 1 2.8.6 


D = D1 + U CR + SO 
n n n n n 


n = N, N-1, 1 


2.8.7 


Note that » Ul^ and that the direct irradiance at each level has been 


Included in D • 
n 
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3. COMPUTATIONAL OVERVIEW 


This section discusses the general flow of the solar radiation code and the 
function of the various subroutines. 

3.1 Initialization 

The input variables for the solar radiation code come from a number of 
sources, depending mainly on their time dependence. This is reflected In the 
schematic flow diagram of fig. 3.1. 

Climatological, or otherwise predetermined, variables enter the radiation 
code by way of DATA statements in the relevant subroutines or may be read in 
from tape at the start of GCM run. 

Slowly varying quantities such as the solar declination angle and the mean 
earth-sun distance are updated at the start of each new Julian day by a call to 
SRDATE. 

Rapidly changing quantities such as temperature, pressure, humidity and 
clouds are calculated at each time step of the GCM and are available to the 
solar radiation code by way of COMMON blocjcs. 

As shown In fig. 3.1 a call to the solar radiation code Involves looping 
over each horizontal grid element, determining the solar zenith angle for that 
element, and calling SOLARl to find the absorbed radiation at each vertical level, 
given a daytime situation. 

3.2 SOLARl 

Subroutine SOLARl controls most of the flow for calculating absorption at a 
single grid point. It also performs the absorption calculations directly In the 
case of a cloud free atmosphere. 

As shown schematically In fig. 3.2, the water vapor amount above each level, 
scaled for pressure and temperature effects. Is determined first. The levels 
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Enter slowly varying input parameters SOMULT, AG and OZONE 
from tape or other permanent storage. 




■j 

i 

.} 

I 


( 


- } 




'1 


i 


\ 

f 

For each new JDAY call SRDATE to update the astronomical 
variables DEC, SIND, COSD, EQT. 

. ___ _ i 

f 

For each time step requiring 
absorption, enter the 
TOFDAY, CL, SHL, TL and PRES 

the calculation of solar 
current values of ! 

S as calculated b/ GCM. 1 

i 




Loop over the entire range of LAT and LON 


Evaluate COSZ 




COSZ -0 => Zero solar absorption 


u. 


COSZ >0 call SOLARl to evaluate the solar 

\ 

absorption by each level of the j 

horizontal grid element. 


J 


Fig, 3.1 Schematic flow diagram of initialization and call to SOLARl 
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Find the scaled water vapor amount above each level. 


Find the cloud optical thickness, if any, at each level. 


If no layers contain cloud. 


ORiG'^' 

OFPco:“;- 


Find the albedo of the surface plus Rayleigh atmosphere outside the 
water vapor absorption bands, using LH(41) and (15). 


Determine the water vapor absorption, neglecting atmospheric 
scattering but including surface reflection, using LH(21). 


Determine surface absorption resulting from radiation transmitted 
through the Rayleigh atmosphere plus that not absorbed by water vapor. 


If cloud is present. 


Call CLOUDY to determine cloud absorption, reflection and transmission. 


Determine surface absorption from the cloud transmission. 


Call 0Z0NE2 to determine ozone amount. 


Determine ozone absorption using LH(8) and (9), allowing 
for reflection from the underlying atmosphere. 


Fin. 3.2 Schemtl'* flow di n i for AR1 








are then examined for cloud content and cloud optical taickneeses assigned baaed 
on cloud type. Should cloud be present at any level the flow turns to subroutine 
CLOUDY for determination of the cloud properties » otherwise absorption due to 
water vapor in Rayleigh atir^sphere is calculated directly. 

SOLARl finally determines the albedo of the earth-troposphere and calculates 
the ozone absorption in the upper levels. 

3.3 CLOUDY 

Subroutine CLOUDY controls the calls to the radiative transfer routines and 
the wavelength integration when clouds are present. 

The effect of clouds in the visible is handled simply by adding all the 
cloud layers together and assuming no atmospheric absorption, so that a single 
call to LAYER and one to ADDER with N=2 to add in the surface albedo is all that 
is needed. 

In the solar infrared, water vapor absorption is included by first deter- 
mining the total intrinsic radiative properties for each layer arising from both 
the cloud and the water vapor, making a call to LAYER and repeating this procedure 
for each layer. ADDER is then called to combine the effects of all the cloud 
layers and the surface, and this process is repeated for each wavelength interval. 


3.4 LAYER 

Subroutine LAYER evaluates the reflection and transmission of a single or 
composite cloud layer, both for unit diffuse irradiance on the cloud top and for 
a given irradiance from a specified solar zenitn angle. 

The cloud properties are scaled using 2.7.9 to 2.7.11 If t > 1. For t < 8, 
the Coakley-Chylek two stream approximation of section 2.7.2 is used- Otlierwise, 
for T ^ 8, the delta Eddington approximation of section 2.7.1 is used, as depicted 
in the schematic fig. 3.4. 

orig::. . ; _ 

OF POOR QL’ALITY 
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f 





Determine absorption of direct radiation 

due to water vapor above clouds. 


5 

t 



Determine the effects of non-conservative scattering in clouds. 



\ 

/ 



Loop over wavelength using the k-distribution. 



1 



Loop over each vertical level. 


i 




Determine the composite cloud-water vapor properties 
and call LAYER. 

j 



i. 


Call ADDER 


to combine all levels plus the surface. 


♦ 


Determine absorption of reflectec 

above the cloi 

radiation by water vapor 
id top. 

\ 

( 

Return to SOLARl 




Fig. 3.3 Schematic flow of subroutine CLOUDY 
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Determine diffuse refU'ttlon and ’•.ransmlsslon using the i 

I Sagan and Pollack two tream. : 

i . I I 


If T >1 Scale the layer properties for a delta 
phase function. 


I * 

I If T ^8 Evaluate U0,D0 for the layer using the j 

j Coakley and Chylek two stream approximation I 


, t 

I " ^ ' t 

If T >8 Evaluate U0,D0 for the layer using the 


I delta Eddington approximation j 


\ 

f 

Return to CLOUDY j 

1 1 

1 ! 


Fig. 3.4 Schematic flow diagram for subroutine LAYER 
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3.5 0Z0NE2. . ADZOME and AWATER 

Subroutine OZONE2 interpolates from a stored data set to find the cliaato- 
logically expected ozone amount above each level, based on the time of year and 
the latitude. 

Subroutine ADDER conblnes the effect of different cloud levels and/or the 
surface. It follows section 2.8 closely. 

Functions AWATER and AOZONE determine the absorption due t? a given path 
of water vapor or ozone, determined from LH(21) and LH(8) , (9) rerreccively . 


3.6 SUBROUTINE SUMM^IES 

3.6.1 AOZONE 

Real Function AOZONE (X) 

This function calculates the fraction of incident radiation that is absorbed 
by an ozone amount of X cm (NTP), using LH(7) and (9). 

There are no external references. 

3.6.2 AWATER 

Real Function AWATER (X) 

This function calculates the fraction of incident radiation that is absorbed 
by X cm of precipltable ^ater vapor, using LH(21) . 


3.6.3 ADDE R 

Subroutine ADDER (NLI) 

COMMON/ ADD/ SOL(IO), UO(IO), DO(IO), RL(IO), TL(IO), UL(IO), DL(IO), ABSL(IO). 
Input Variables 

NLI - Integer - The number of layers to be combined. 

2 < NLI < 10 

SOL - Real Array - Vertical component of the direct solar beam crossing 


each level. 




UO ~ Real Array - Upwelllng Irradiance at the top of each layer due 

to thi- direct solar beam on that layer • 

DO > Real Array - Downwelling Irradiance at the base of each layer. 

due to the direct solar beam. 

RL - Real Array - Reflectivity of each layer to diffuse irradiance. 

TL - Real Array - Transmissivity of each laver to diffuse irradiance. 

Output Variables 

UL - Real Array - The total upwelllng irradiance through 

each level. 

DL - Real Array - The total downwelling irradiance thr< ugh 

each level. 

ABSL - Real Array - The solar radiation absorbed by each layer. 

Notes 

1. This subroutine is called from CLOUDY for each spectral interval to 
evaluate the absorption due to interactions between layers. 

2. The computer variables generally have names which compare uirectiy 
with the symbols of section 2.8. 

3. The first part o^ ADDER is for NL1=2, which Cakes a particularly simple 
form. The section is entered to evaluate the total cloud effect In the 
nonabsorbing portion of the spectrum. 


3.6.4 CLOUD Y 

Subroutine CLOUDY (RCLOUD, NLAY, NTOP) 
CONMON/RADCOIl/ 
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COMMON/CLDCCW/ 
CCMMON/ADD / 


Input Variables 


NIAY 

- Integer 

NTOP 

- Integer 

TAUL(16) 

- Real Array 

FSCAT 

- Real 

COSZ 

- Real 

RSURF 

- Real 

AG 

- Real 

SUALE(16) 

- Real Array 

COSMAG 

- Real 

SCOSZ 

- Real 

SWIl'16) 

- Real Array 


Number of layers in model. 

Layer number of the highest layer containing cloud 
Cloud optical thickness of each layer . (CLDCOM) 

Fraction of incident radiation (visible) that is 
not subjected to atmospheric absorption (CLDCOM) . 

Coslae of the solar 7.enith angle (RADCOM) • 

Surface albedo with respect to direct beam (CLDCOM) . 
Surface albedo with respect to diffuse radiation(CLDCOM) 
Scaled water vapor amount above each level (CLDCOM) . 
Refraction magnification factor (CLDCOM) • 

Solar irradiance normal to top of atmosphere (CLDCOM) . 
Scaled water content within each layer (CLDCCW) . 


Output Variables 

TOPABS - Real 

AL(16) - Real Array 

RCLOUD - Real 


Absorption above level 1 (CLDCOM). 
Absorption in each layer (CLDCOM) . 
Cloud albedo in the visible. 


Notes 

1. This subroutine first evaluates the effect of all the cloud layers in the 
visible by a single call to LAYER and one call to ADDER . 

2. Vhe variables required for ALDER a transferred using C(M!ON/ADD/. 

3* This subroutine is called by SOLARl and in turn makes use of LAYER and ADDER. 



3.5.5 LAYER 




Subroutine LAYER (SCOSZ» COSZ, TAU» PIO, SOl» UP, DN, RL, TL) 

Input Arguments 

SCOSZ - Real - Vertical component of direct solar Irradiance on top 
surface of layer. 

COSZ Real * cosine of the solar zenith angle. 

TAU - Real ~ Layer optical thickness. 

PIO * Real * Single scatter albedo. 

Output Arguments 

SOI - Real - Vertical component of transmitted direct solar irradiance. 
UP - Real - Irradiance reflected from the direct beam by the layer. 

DN - Real - Irradiance transmitted through the layer from the direct 

beam. 

RL - Real - The reflectivity of the layer to diffusely incident 
radiation. 

TL - Real - The transmissivity of the layer to diffusely incident 
radiation. 

Notes 

1. This subroutine is called by CLOUDY to evaluate the reflection and 
transmission characteristic of any layer, such as a cloudy layer, that 
contains scatterers. 

2. The output from LAYER is utilized by ADDER. 

3. Subroutine LAYER is self-^contalned, with no external references. 

4. The computer variables generally have names which compare directly 
with the symbols of section 2.7. 
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3.6.6 0Z0NE2 

Subroutine 0Z0NE2 (NUYOl, XDAY, XLAT) 
COMfON/RADCOM/ 

COMfON/CLDCOl/ 


Input Variables 

NLAYOl - Integer - The nund>er of pressure levels above which ozone 

is to be considered. SOLARl sets NLAYOl > 6. 
XDAY - Real - CARP reference day. XDAY • JDAY + 63. 

XLAT - Real - Ladltude in degrees. 

PLE (16) - Real Array - in RADCOM. 

- the pressure (mb) at the top of each level of 

the GCM. 


Output Variables 

OZALE (16) - Real Array - in CLD(X)M. 

- the amount of ozone above the top of each model 
layer. 


Notes 

1. This subroutine is called once per grid point via SOLARl to calculate 
the ozone amounts based on climatology above each level of the model. 

2. There are no external references. 

3. The subroutine first finds the total ozone (TOTOCM) by interpolating 
from quarterly climatological values based on latitude. It then uses 
TOTOCM to set up a fine scale vertical distribution of ozone from which 
it interpolates to obtain the values of OZALE. 
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3.6.7 SQLARl 


Subroutine SOLARl (NLAY, XIjAY, XLAT) 
C(MfON/RADC(»l/ 

C(MiON/CLDC(»l/ 


Input Variables 


NLAY 

- Integer 

XJJAY 

- Real 

XLAT 

- Real 

COSZ 

- Real 

RSURF 

- Real 


- The number of atmospheric layer:;. 

- GARP reference day. XDAY - JDAY + 63. 

- Latitude in degrees. 

- Cosine of solar zenith angle. (RADCCW) 

- Surface albedo for direct radiation from 


COSZ. (CLDCOM) 


AG - Real - conventional surface albedo. 

so -- Real - Solar irradiance at top of atmosph re (RADCOM) . 

PLE (16) - Real Array - pressure (mb) at ^ach atmospheric level 

(RADCOM) . 


TL 

(i:» 

- Real Array 

SUL 

(15) 

- Real Array 

CL 

(15) 

- Real Array 


Temperature at mid-pt. of each layer (RADCOM). 
Humidity of each layer ( RAT COM) . 
cloud information, coded to indicate presence 
and type. (RADCOM) 


Output Variables 


AS (15) -- Real Array 

SG - Real 

SWALE (16) - Real Array 
SWIL (16) - Real Array 


absorption in each atmospheric layer (RADCOM) 
absorption at surface. (RADCOM) 
scaled water vapor above each level. (CLDCOM) 
scaled water vapor amount in each layer. (CLDCOM) 
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Notes 


1. The voter vapor la scaled for tei^>erature and pressure effects as 
In LH(2A). 

2. The cloud optical thickness In each level is found by decoding the 
array CL and follovs section 2.6. 

3. The cloud free situation is treated directly, fli stly for the non- 
absorbing Rayleigh atoosphete, nd secondly for water vapor absorp- 
tion in the absence of scattering. 

4. This subroutine is called at each grid point to evaluate solar 
absorption. It calls in turn CLOUDY, if clouds a.*e pres'^at, and 
0Z0NE2 for ozone amounts. 

3.6.8 SRDATE 

Subroutine SRDATE (JDAY, S(MILT, SO, DEC, SIND, COSD, EQT, ALBEDO) 


Input Arguments 

JDAY 

- Integer 

• The Julian day nunA>er» * 1 on Jan. 1st. 

SOMULT 

- Real 

- A multiplier of the solar constant. 
SOMULT » 1 for units of ly/day. 

SORJLT - 0.4844 for units of W/M^ . 

SOMULT can also be varied to perturb the 
solar constant. 

Output Arguments 

SO 

- Real 

- The extraterrestrial solar irradiance. 

DEC 

- Real 

- The solar declination angle (radians). 

SIND 

- Real 

- SIN(DEC). 

COSD 

- Real { 

1 

- COS(DEC). 

EQT 

- Real 

- Equation of tine (radians). 

ALBEDO (46,72) 

- Real Array 

- Conventional surface albedos. 



Notes 


1. This subroutine is called once per day to update slowly varying 
astronomical variables and surface albedos. 

2. A solar constant of 1368 W/m^ is used. 

3. The equations used are discussed in section 2.1. 

4. The surface albedos are read in from unit 10 as lALB (46,72) - Integer 
*2 in free fomat. 

Calling Routine: Main Program 

External References: Logical Unit 10. 
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4. Comparative Results 


Results from the new code for the transmittance of a single cloud layer as 
a function of sun angle are compared In Fig. 4.1 (taken from Davies* 1980) with 
results from the Lads and Hansen (LH) code and precise results from a doubling 
model (Hansen* 1971). While the new code compares well with the doubling 
model* the LH version underestimates transmission for high solar elevations 
and overestimates it for low solar elevations. 

In addition to specific case comparisons between the old and new codes* 
which generally show large differences* it is of interest In the GCM context 
to compare averaged results to see if the differences might be significantly 
reduced. The GCM output was therefore analyzed off-line to provide such 
averages. The input parameters to the radiation routines obtained during a 
given GCM run were stored and the absorbed radiation was recalculated with 
both the old and new codes. Since the new code did not interact with other 
aspects of the GCM, a direct comparison of code differences was therefore 
possible* circumventing the immediate effects of model noise. 

The off-line recalculation oil absorbed solar radiation using the LH 
code was of course redundant* the values agreeing with chose used in the 
original on-line calculations of the complete GCM. The recalculation did* 
however* serve as a diagnostic check and conveniently reformatted the old 
results . 

Figs. 4. 2-4.4 show results from the new code (solid lines) and the 
differences with the LH code (dashed lines: LH result-new result). The 

absorbed solar radiation has been zonally averaged and also averaged for the 
mouth of February* for a GCM run which was Initialized on January 1. Fig. 4.2 
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shows the radiation absorbed by both the surface and ataosphere. Fig. 4.3 and 
Fig. 4.4 show, respectively, the surface and atmospheric absorption. The LH 
code is seen to systematically overestimate absorption at almost all latitudes. 
The differences vary sooiewhat with latitude, depending on sun angle and cloud 
cover. On a globally averaged basis, the overestimate Is 7Z in total absorption, 
4% In surface absorption and 17% In atmospheric absorption. The corresponding 
results for zonally averaged planetary albedo are shown In Fig. 4.5. The LH 
code systematically underestimates albedo with the differences Increasing at 
high latitudes. Fig. 4.6 shows similar results, but for an average over April. 
If the time average Is relaxed and zonal averages at a single time step are 
compared, as In Fig. 4.7 for mid-February, the functional dependence on latitude 
is less regular, but the systemaf'.c differences persist. 

Finally, since the largest percentage differences were in atmospheric 
absorption, the vertical distribution of absorbed solar radiation by model 
level Is presented In Fig. 4.8 for the global average and a zonal average at 
46“S. The two codes show their greatest differences in levels 4 and 5, pre- 
sumably due to the effects of thick cloud. 
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Fig. 4 


2 Zonally averaged solar radiation absorbed by the surface plus 
atmosphere. February average. 

New code results. — -- — 


LH code - New code. 
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Fig. 4.4 Sane as Fig. 4.2 but for acmospherlc absorption only. 
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Pig. 4.5 Zonally av«r«g«d planetary albedo. February average. 
Ne«r code. — LH code. 




35 



PLANETARY ALBEDO 
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Fig. 4.7 S«ae as Fig. 4.5 but for single tine step at mid-February. 
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Fig, 4.8 Atmospheric absorption by cnodel level* February averages. 
New code* — — — I.H code. 

(a) Global average (b) Zonal average <»t 46 ® S 
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5. S u—ary 


The new solar radiation code described in this report yields significantly 
different values of absorbed solar radiation by the atmosphere and surface when 
compared to the original Lacis and Hansen code* These differences appear 
greatest for conditions of large solar zenith angle, and consequently, depend 
on latitude and season. 

While the new code follows the physics of radiative transfer quite faith- 
fully, as evidenced by comparison with even more precise (yet slower) models, 
there undoubtedly remains some room for improvement. A spectrally dependent 
surface albedo, and the effect of a background aerosol, should be included, but 
the most substantial improvement will be in the treatment of cloud propoerties, 
especially the treatment of fractional cloud cover* A substantial amount of 
original research on the radiative properties of multi-layered broken clouds Is 
first needed, however, before they can be successfully considered in a GCM. 
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REAL FUNCTION AOZONE(X) 

C 

€*★* THIS IS THE UCIS-HANSEN PARAMETERIZATION OF OZONE ABSORPTION. 

C 

IF(X.LE.O. )GOTO 1 

AOZONE-.O2118O*X/(1.0+0.OA2*X+0.O0O323*X*X)+l.O8i73*X/(l.OO+ 

1 138 . 57*X) **0 . 805+0 . 0658*X/ (1 . 00+ (103 . 63*X) **3) 

RETURN 

C 

1 AOZONE-0.0 
C 

RETURN 

END 

REAL FUNCTION AWATER(X) 

C 

C*** THIS IS THE LACIS-HANSEN PARAMETERIZATION FOR WATER VAPOR 
ABSORPTION. 

1F(X.LT.0. )GOTO 1 

AWATER-2 . 9*X/ ( ( 1 . + 1 41 . 5 1 *X) ** . 635+5 . 925*X^ 

RETURN 
C 

1 AWATER-0.0 

C 

RETURN 

END 

SUBROUTINE ADDER (NLl) 

C 

C COMBINES THE SOLUTIONS OBTAINED PREVIOUSLY FOR ISOLATED LAYERS, 
C BY SUCCESSIVE APPLICATION OF THE INTERACTION PRINCIPLE 
C 

C DESIGNED BY ROGER DAVIES (PURDUE UNIVERSITY), 1981. 

C 

C NL IS THE NUMBER OF ISOLATED LAYERS TO BE ADDED, INCLUDING 
C THE SURFACE, AND NLl IS NL+1. 

C 

C ASSUMES THAT NLl IS BETWEEN 2 AND 10, AND THAT U0,b0,S0L,RL 
C AND TL HAVE BEEN SET UP BY CLOUDY. UL.DL AND ABSL ARE FOUND 

C HERE. 

C 

C THE TEN'S BELOW ARE NLAY+1, AND SHOULD BE CHANGED IF THE 
C NUMBER OF MODEL LAYERS IS CHANGED. 

C 

REAL M(10) ,CR(10) ,VM(10) ,VP(10) 

C 

COMMON/ADD/SOL (10) ,U0(10) ,D0(10) ,RL(10) ,TL(10) , 

1 UL(10),DL(10),ABSL(10) 

C 

NL-NLl-1 

M(2)-1./(1.-RL(1)*RL(2)) 

IF(NL.GT.l )GOTO 10 

UL (2) - (UO (2) +D0 (2) *RL (2) ) *M (2) 

UL(1)-U0(1)+UL(2)*TL(1) 

DL(2)-S0L(2)+D0(2)+UL(2) *RL(1) 

ABSL(1)-S0L(1)+UL(2)-UL(1)-DL(2) 

ABSL(2)-DL(2)-UL(2) 

RETURN 

C 

10 CR(l)-0.0 


SWDOOOlO 

SWD00020 

SWD00030 

SWD00040 

SWD00050 

SWD0006O 

SWD00070 

SWD00080 

SWD00090 

SWDOOlOO 

SWDOOllC 

SWD00120 

SWD00130 

SWD00140 

SWD0015C 

SWD00160 

SWD00170 

SWD00180 

SWD00190 

SWD00200 

SWD00210 

SW00022C 

SWD00230 

SWD00240 

SWD00250 

SWD00260 

SWD00270 

SWD00280 

SWD00290 

SWD0030C 

SWD00310 

SWD00320 

SWD00330 

SWD00340 

SWD00350 

SWD00360 

SWD00370 

SWD00380 

SWD00390 

SWD00400 

SWD00410 

SWD00420 

SWD00430 

SWD00440 

SWD00450 

SWD00460 

SWD00470 

SWD00480 

SWD00490 

SWD00500 

SWD00510 

SWD00520 

SWD00530 

SWD00540 

SWD00550 

SWD00560 

SWD00570 

SWD00580 

SWD00590 
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3 

C 


r 

C*** 


CR(2)-RL(1) 

DO 1 J-2.nl, 

CR (J+D-RL f J) -^TL( J) *TL(J) *CR ( J) *M ( J) 

::(j-M)-i./(i.-CR(j+i)*RL(j+i)) 

VP (2) -DO (2) 

VP(l)-0.0 
VH(1)-U0(1) 

DO 2 J-l.NL 

VP ( J-*-!) -TL ( J) *M ( J) * (VP ( J) +CR ( J) *UO( J) ) +D0 ( J+ 1) 
VM(J+1)-M(J+1)*(VP(J+1)*RL(J+1)+U0(J+D) 

UI.(NLl)-VIt(NLl) 

DL (NLl) -VP (NLl) +UL (NLl) *CR (NLl) +SOL (NLl) 

ABSL (NLl ) -DL (NLl) -UL (NLl ) 

DO 3 J-l.NL 
K-NLl-J 

UL (K) - VM (K) +UL (K+ 1 ) *TL (K) *M (K) 

DL (K) -VP (K) +UL (K) *CR (K) +S0L (K) 

ABSL (K) -DL (K) -UL (K) - (DL (K+ 1 ) -UL (K+ 1 ) ) 

RETURN 
END 

SUBROUTINE CLOUDY (RCLOUD. KLAY , NTOP) 

THIS SUBROUTINE FIRST CREATED BY R. DAVIES ll/U/79 

IT EVALUATES THE ABSORPTION DUE TO WATER VAPOR IN A CLOUDY 
ATHOSPHEHE. CALLING LAYER AND ADDER IN THE PROCESS. 

THE EXISTENCE OF CLOUD IN AT LEAST ONE LEVEL BELOW THE TOP 
LEVEL IS ASSUMED. 

SWD00920 

COMMON/RADCOM/AS(15) ,RE(16) ,PL(15) .PLE(16) ,PLK(15) ,PLK£C16) .TT(15)SWD00930 
*,TLE(16) ,TG.TH(15) ,SHL(15) .SHLE(16) .SHG,CLOUD(15) .COSZ.SO.SC.CXL, SWD009A0 
* OCEAN. ICE. SNOW SWD00950 

SWD00960 

LOGICAL OCEAN, ICE, SNOW SWD00970 

SWD00980 

C0MM0N/CLDC0M/SWALE(16) ,SWIL(15) ,AL(16) ,TAUL(16) ,0ZALE(16) ,TOPABS,SWD00990 
1 C0SMAG,SC0SZ,FSCAT,RSURF,AG SWDOIOOO 

SWDOlOlO 
SWD01020 
SWD01030 
SWUC1040 
SWD01050 
SWD01060 
SWD01070 
SW001080 
SWD01090 
SWDOllOO 

C SWDOlllO 

C EVALUATE THE CONSERVATIVE PORTION OF THE CLOUD ALBEDO. SWD01120 

SWD01130 
5WD011A0 
SWD01150 
SWDOllbO 
SWD01170 
SWDOllSO 

CALL LAYER (SOL(l) ,COSZ.TAU, 0.99 '99,501(2) ,U0(1) ,D0(2; ,RLU) T1 ; '))SWD01 ’QO 


SWD00600 

SWD00610 

SWD00620 

SWD00630 

SWD00640 

SWD00650 

SWD00660 

SWD00670 

SWD03680 

SWD00690 

SWD00700 

SWD00710 

SWD00720 

SWD00V30 

SWD00740 

SWD00750 

SWD00760 

SWD00770 

SWD00780 

SWD00790 

SWD00800 

SWD00810 

SWD00820 

SWD00830 

SWD00840 

SWD008S0 

SWD00860 

SWD00870 

SWD00880 

SWD00890 

SWD00900 

SWD00910 


COMMON/ADD/SOL (10) ,U0(10) ,DO(10) ,RL(10) ,TL(10) , 

1 UL(10).DL(10),ABSL(1D) 

DIMENSION FK(5).XX(5) 

DATA FK/0. 107, 0.104,0.073,0.044,0.025/ 

DATA XK/0. 005. 0.041. 0.416,4.752, 72. 459/ 

DATA NFK/5/.EPS/1./ 

NLAYl-NLAY+1 

NCLEAR“NT0P-1 

EVALUATE THE CONSERVATIVE PORTION OF THE CLOUD ALBEDO. 
TAU-O.C 

DO 1 N-NTOP,NLAY 
TAU-TAU+TAUL(N) 

S0L(1)-FSCAT 


44 



t 


OMGIMAL' PAGE IS 


c 

OF POOR QUALITY 

SWD01200 


i* 


U0(2)-RSURF*S0L(2) 

SWD01210 


1 

RL(2)-AG 

SWD01220 


t 

C 


SWD01230 



CALL AIH)ER(2) 

SWD01240 



c 


SWPO1250 




RCLOUD-UL(l)/FSCAT 

SWiX)1260 



c 


SWD01270 



r 

V 

EVALUATE ABSORBED INCIDENT FLUX ABOVE CLOUDS. 

SW001280 



c 


SWD01290 



V#-SWALE(l)’^COSMAG 

SWD01300 


. 1 - 
t 

DB-AWATER (W) 

SWD01310 


f 

T0PABS«DB*SC0S2 

SWD01320 


c 


SWD01330 



DO 2 N-l.NCLEAR 

SWD01340 


t 

W-SWALE(N+l)*COSMAG 

SWD01350 


f 

1 

DA-AWATER (W) 

SWD0I360 



AL(N)-(DA-DB)*SC0S2 

SWD01370 


' 2 

DB-DA 

SWD01380 


C 


SWD01390 


C 

EVALUATE ABSORPTION WITHIN AND BELOW CLOUDS. 

SWD01400 


C 


SWD01410 



DO 3 N-NTOP.NUYl 

SWD01420 


: 3 

AL(N)-0.0 

SWD01430 


c 


SWD01440 


1 c 

LOOP OVER WAVELENGTH USING THE K-DISTRIBUTION. 

SWD01450 


i ^ 


SWD01460 


I 

DO 4 K-l.NFK 

SWD01470 



WK«W*XK(K) 

SWD01480 



WK-AMIN1(WK,75.0) 

SWD0I490 


I 

SOLd)-EXP(-WK) *FK(K) *SC0S2 

SWD01500 


r 

IF(S0L(1).LT.EPS )rOTO 4 

SWD01510 


1 

DO 5 N-NTOP.NLAY 

SWD01520 



M-N-NCLEAR 

SWD01530 



TAUAB-XK (K) *SWI L (N) 

SWD01540 


t 

i 

TAU-TAUL(N)+TAUAB 

SWD01550 


( 

PI0-TAUL(N)/TAU 

SWD01560 


( 

IF(PI0.GT.0.1 )GOTO 7 

SWD01570 



ARG-AMINl (TAUAB*COSHAG, 75 . 0) 

SWD01580 

1 


SOL (M+ 1) -SOL (H) *EXP (-ARG) 

SWD01590 

1 


U0(M)-0.0 

SWD01600 

1 

1 

f 

DC(M+l)-0.0 

SWD01610 


RL(M)-0.0 

SWD01620 

1 


ARG-A.'lINl (1 . 66*TAUAB ,75.0) 

SWD01630 

1 

! 

TL(H)-EXP(-ARG) 

SWD01640 

1 

i 

GOTO 5 

SWD01650 

1 

, f c 


SWD01660 

I 

! ! . ^ 

CALL LAYER (SOL (M) , C0S2 , TAU , PI 0 , SOL (H+ 1 ) , UO (M) , DO (M+ 1 ) 

,RL(M),SWD01670 

1 

i 

1 TL(H)) 

SWD01680 


1 ^ 

CONTINUE 

SW001690 


\ ! 

M-NUYl-NCLEAR 

SWD01700 

1 

! 


U0(M)-RSURF*S0L(M) 

SWD01710 

i 

a 

I 

RL(M)-AC 

SWD01720 

* 

1 

c 


SWD01730 




CALL ADDER (N) 

SWD01740 



C 


SWD01750 




DO 8 N-NTOP.NUYl 

SWD01760 



M-N-NCLEAR 

SWD'J1770 


8 

AL(N)-AL(N)+ABSL(M) 

SWD01780 


C 


SWD01790 

i 


I 



n n 


I 


C EVALUATE ABSORPTION OP REFLECTED FLUX ABOVE THE CLOUDS. 

C 

FKK-UL(l) 

N-NTOP 

9 IF(PKK.LT.EPS )GOTO A 

N-N-1 

ARG-AHIN1(1.66*SWIL(N)*XK(K) ,75.0) 

DA- ( 1 . 0-EXP (-ARC) ) *FKK 
AL(N)-AL(N)+DA 
IF(N.GT.l )GOTO 9 

ARG-AMINI (1 . 66*SHIL (1) *XK (K) ,75.0) 

TOPABS-TOPABS+FKK* (1 . 0-EXP (-ARC) ) 

A CONTINUE 

C 

RETURN 

END 

SUBROUTINE UYER (SCOSZ , COSZ , TAU , P 10 , SOI , UP , DN , RL , TL) 

C 

C*** THIS CALCULATES THE REFLECTION, ABSORPTION, AND TRANSMISSION 
C OF AN ISOLATED CLOUD LAYER, USING THE DELTA- EDDINGTON OR 
C TWO-STREAM APPROXIMATIONS, DEPENDING ON LAYER THICKNESS. 

C 

C SCOSZ IS THE COSINE-WEIGHTED INCIDENT SOLAR FLUX FOR THE LAYER. 
C COSZ IS THE COSINE OF THE SOLAR ZENITH ANGLE. 

C TAU IS THE LAYER’S OPTICAL THICKNESS. 

C PIO IS THE SINGLE-SCATTERING ALBEDO FOR THE LAYER. 

C 

C SOI IS THE DIRECT FLUX TRANSMITTED BY THE LAYER. 

C UP IS THE UPWARD DIFFUSE IRRADIANCE AT THE TOP OF THE LAYER, 

C IN UNITS OF SCOSZ. 

C DN IS THE DOWNWARD DIFFUSE IRRADIANCE AT THE BASE OF THE LAYER, 
C IN UNITS OF SCOSZ. 

C RL IS THE REFLECTANCE OF THE LAYER. 

C TL IS THE DIFFUSE TRANSMITTANCE OF THE LAYER. 

C 

REAL KP,KS,IOO,IOZ,LAMDA,LAMDA2 
C 

REAL BA(17' ,BB(12) 

C 

DATA PI/3. 1A159265A/.F/0. 7225/ 

C 

DATA BA/0.5, .32, .2, . lA, . 11 , .09, .07, .06, .05, .0A5, .04, .04/ 

C 

DATA BB/. 5, .47, .43, .38, .36, .32, . 29, . 265 , . 24, . 215 , . 19, . 19/ 

C 

DATA G/0.85/ 

C 

SECZ-l./COSZ 

ZO-TAU 

A-PIO 

C 

C FIRST CALCULATE THE DIFFUSE REFLECTANCE AND TRANSMITTANCE, 

C USING THE SAGAN-POLLACK TWO-STREAM APPROXIMATION. 

A-1 IMPLIES CONSERVATIVE SCATTERING (NO ABSORPTION). 

IF(A.LT.1.0)G0TO 12 

TL-l./(l.-HO.13*Z0) Uk.G... cJ 

RL-I.-TL of poor quality 

c 


SWD01800 

SWD01810 

SWD01820 

SWD01830 

SWD01840 

SWD01850 

SWO01860 

SWD01870 

SWD01880 

SWD01890 

SWD01900 

SWD01910 

SWD01920 

SWD01930 

SWD01940 

SWD01950 

SWD01960 

SWD01970 

SWD01980 

SWD01990 

SWD02000 

SWD02010 

SWD02020 

SWD02030 

SWD02040 

SWD02050 

SWD02060 

SWD02070 

SWD02080 

SWD02090 

SWD02100 

SWD02110 

SWD02120 

SWD02130 

SWD02140 

SWD02150 

SWD02160 

SWD02170 

SWD02180 

SWD02190 

SWDQ2200 

SWD02210 

SWD02220 

SWD02230 

SWD02240 

SWD02250 

SWD02260 

SWD02270 

SWD02280 

SWD02290 

SWD02300 

SWD02310 

SWD02320 

SWD02330 

SWD02340 

SWD02350 

SWD07360 

SWD02370 

£.^007380 

C .•'■"'90 
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n n 


t .. 




ORIGINAL PAGE (S 

SCATTERING PLUS ABSORPTION. OF POOR QUALITY 

12 U2-(I.-A*G)/(l.-A) 

U-SQRT(U2) 

DL-SQRT (3 . * ( 1 . -A) * ( 1 - -A*G) ) 

DS-DL*ZO 

IF(DS.GT.75.)DS-75. 

DO-EXP (-DS) 

D02-D0*D0 

DEN- (U+1.)*(U+1.)-(U-1.)*(U-1.)*D02 
RL- (U2- I . ) * ( 1 . -D02) /DEN 
TL-4. *U*Do/DEN 
C 

C NOW CALCUUTE the reflection and transmission DUE TO THE 
C DIRECT BEAM. 

C 

C IF TAU IS LESS OR EQUAL TO 1, WE USE THE COAKLEY-CHYLEK 
C TWO-STREAM APPROXIMATION, WITHOUT SCALING. IP TAU IS 
C GREATER THAN 1 AND LESS THAN OR EQUAL TO 8. USE THE 

C COAKLEY-CHYLEK TWO-STREAM APPROXIMATION, WITH SCALING. 

C FOR TAU GREATER THAN 8. USE THE DELTA-EDDINGTON 

C APPROXIMATION. 

C 

13 IF(TAU.GT.l. )GOTO 1 
C 

P-SECZ 
GOTO 2 
C 

1 Z0-TAU*(1.-P10*F) 

G-O.A5945946 
A-PI0*0. 2775/ (1 . -PIO*F) 

2 DS-ZO*SECZ 
1F(DS.GT.75.)DS-75. 

EO-EXP(-DS) 

S01-SCOSZ*E0 

c 

IF(TAU.GT.8. )GOTO 3 
C 

BETA-0. 5-0. 4375*G 
IR-10*COSZ 
C 

IF(TAU.GT.l. )G0T0 4 
C 

BETA2-BA ( I R+ 1 ) - ( 1 0 . *COSZ- IR) * (BA ( IR+ 1 ) -BA ( lR+2) ) 
BEiAl-BETA2 
GOTO 5 
C 

4 BETA2-BB (IR+1) - (10. *COSZ-IR) * (BB (IR^l) -BB (IRi^2) ) 
BETAl-BETA 

C 

P-2.0 

C 

5 LAMDA2-P*P*(1.-A)*(1.-A+2.*A*BETA1) 

EP— A* (P*A*BETA1 ♦ ( 1 . -BETA2) * (P* ( 1 . - A) <^SECZ) ) 

EM— A* (P*A*BETA1+BETA2* (P* (1 . -A) -SECZ) ) 

T-SECZ*SECZ-LAMDA2 

GP-EP/T 

GM-EM/T 

c 

C CHECK FOR CONSERVATIVE SCATTERING. 


SWD02400 

SWD02410 

SWD02420 

SWD02430 

SKD02440 

SWD02450 

SWD02460 

SWD02470 

SWD02480 

SWD02490 

SWD02500 

SWD02510 

SWD02520 

SWD02530 

SWD02540 

SWD025SU 

SW002560 

SWD02570 

SWD02580 

SWD0259C 

SWD02600 

SWD02610 

SW002620 

SWD02630 

SWD02640 

SWD02650 

SWD02660 

SWD02670 

SWD02680 

SWD02690 

SWD02700 

SWD02710 

SWD02720 

SWD02730 

SWD02740 

SWD027S0 

SWD02760 

SWD02770 

SWD02780 

SWD02790 

SWD02800 

SWD02810 

SWD02820 

SWD02830 

SWD02840 

SWD02850 

SWD02860 

SWD02870 

SWD02880 

SWD02890 

SWD02900 

SWD02910 

SWD02920 

SWD02930 

SWD02940 

SWD02950 

SWD02960 

SWD02970 

SWD02980 

SWD02990 




onooujn 


I F ( (1 . -A) . LE . 0 . 0) GOTO 6 OF POOR QUALITY 

LArtDA-SQRT(LAMDA2) 

DS-LAMDA*Z0 

IF(DS.GT.75.)DS-?5. 

DO- EXP (-DS) 

Al-l.-A+A*BETAl 

Bl-LAMDA/P 

DEN-Al+Bl- (Al-Bl) *D0*D0 

AP- (GP* (Al-Bl) *D0-GM*A*BETA1*E0) /DEN 

AM-AP* (Al+Bl) / (A*BETA1) 

BP— GP-AP*D0 

BM-BP* (Al-Bl) / (A*BETA1) 

UP- (AM*DO+BH+CM) *SCOSZ*SECZ 
DN- (AP+BP^D0+GP*E0) *SC0SZ*SECZ 
GOTO 10 
C 

6 SP— GP 

SM- (-GP*P*BETAl *Z0-GM*E0) / (1 . +P*BETAl *Z0) 

rho-p*betai*(sm-sp) 

UP- (GM+SM) *SCOSZ*SECZ 
DN- (GP*EO+RHO*ZO+SP) ’*SCOSZ*SECZ 
GOTO 10 

T-l.-A 

CHECK FOR CONSERVATIVE SCATTERING. 

IF(T.LT.1.E-10 )GOTO 8 

r}L2-3.*(l.-A*G)*T 
DL-SQRT(DL2) 

DS-DL*Z0 

IF(DS.GT.75.)DS-75. 

DO-EXP (-DS) 

TEMP-SECZ*SECZ-DL2 
FSUN-3 . *COSZ*A*G/ (4. *PI) 

AL3— 3 . *A* (1 . +G*T) / (4 . *PI*TEMP) 

H-1.5*(1.-A*G) 

Ul-H-DL 
Vl-H+DL 

Wl— (H+SEC2) *AL3-FSUN 
W2— EO* ( (H-SECZ) *AL3-FSUN) 

DET- V 1 * V 1 -U 1 *U I *D0*D0 
ALl- (V1*W2-U1*W1*DG)/DET 
AL2-(V1*W1-U1*W2*D0)/DET 
100-AL1*D0+AL2+AL3 
I OZ-ALl +AL2*DO+AL3*EO 
UP-2 . *P1 *I00*SC0SZ*SEC2 
DN-2.*PI*I0Z*SC0SZ*SECZ 
GOTO 10 

AL-COSZ* ( (3 . *C0SZ -2 . ) *E0-3 . ’"COSZ-2 . ) / (4 . *P I ) 
AL-AL/(Z0+4./(3.*(l.-G))) 

BETA-COSZ* (3 . *COSZ+2 . ) / (4 . *P ] ) +AL/ ( 1 . 5* ( 1 . -G) ) 

UP— 1 . 5*C0SZ*C0SZ-* 2 . *P I *BETA 
DN— 1 , 5*COSZ*COSZ*EO+AL*ZO*2 . *P1 +2 . *P1 *BETA 

SAFETY NET. DON'T ALLOW NEGATIVE RADIATION. 


SWD03000 

SWO03010 

SWD03020 

SWD03030 

SWD03040 

SHD03050 

SWD03060 

SWD03070 

SWD03080 

SHD03090 

SWD03100 

SWD03110 

SWD03120 

SWD03130 

SWD03140 

SWD03150 

.SWD03160 

SWD03170 

SWD03180 

SWD03190 

SWD03200 

SWD03210 

SWD03220 

SWD03230 

SWD03240 

SWD03250 

SWD03260 

SWD03270 

SWD03280 

SWD03290 

SWD03300 

SWD03310 

SWD03320 

SWD03330 

SWD03340 

SWD03350 

SW003360 

SWD03370 

SWD03380 

SWD03390 

SWD03400 

SWD03410 

SWD03420 

SWD03430 

SWD03440 

SWD03450 

SWD03460 

SWD03470 

SWD03480 

SWD03490 

SWD03500 

SWD03510 

SWD03520 

SWD03530 

SWD03540 

SWD03550 

SWD03560 

SWDC3570 

SWD03580 

590 
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c 

10 


IF(UP.LT.O.) UP-0.0 
1F(DN.LT.0.) WO.O 
IF(RL.LT.O.) RL-0.0 
IFdL.LT.O.) TL-0.0 


0R1G!MAL PA'= J IS 
OF POOft QUALITY 


C 

C 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


c 

c 

c 

c 


SHD03600 
SWD03610 
SWD03620 
SWD03630 
SW003640 
SMO036S0 
SWO03660 
SWD03670 
SHD03680 
SWD03690 
SMD03700 
SWD03710 
SWD03720 
SWD03730 
SWD03740 

C0MM0N/RADC0M/AS(15) ,RE(16) ,PL(15) ,PLE(16) .PLK(15) ,PLKE(16) .TL(15)SWD03750 
*,TLE(16) ,TG,TH(15) .SHLOS) ,SHLE(16) ,SHG.CL0UD(15) .COSZ.SO.SG.CXL, SWD03760 
* OCEAN, ICE, SNOW SWD03770 

SWD03780 

C0MM0N/CLDC0M/SWALE(16) ,SWIL(15) ,AL(16) ,TAUL(16) ,02ALE(16) ,TOPABS.SWD03790 

SWD03800 
SWD03810 


RETURN 

END 

SUBROUTINE 0Z0NE2 (NUYOl ,XDAY ,XUT) 

C 

C*** LACIS-HANSEN OZONE ROUTINE. USES CLIMATOLOGICAL OZONE 
C DISTRIBUTION, INTERPOLATED BY LATITUDE, SEASON, AND ALTITUDE. 


LOGICAL ICE, OCEAN. SNOW 


1 COSKAG,SCOSZ,FSCAT,RSURF,AG 

REAL OLJAN(19) ,OLAPR(19) ,OLJUL(19) ,0L0CT(19) 

REAL OCM22(23) ,OCM30(23) ,0CM38(23) ,OCMA6(23) ,PROCM(23) 
REAL T0T0Z(4) ,CDATE(6) ,OCMXX(23) 


SWD03S20 

SWD03830 

SWD03840 

SWD03850 

SWD03860 

SWD03870 

SWD03880 

SWD03890 

SWD03900 


TOTAL OZONE AMOUNTS AS A FUNCTION OF LATITUDE AND SEASON. 

VALUES ARE GIVEN EVERY 5 DEGREES OF LATITUDE FROM THE EQUATOR 
TO 90 DEGREES. SYMMETRY ACROSS THE EQUATOR IS ASSUMED. 

OLJAN IS FOR NORTHERN HEMISPHERE JANUARY AND SOUTHERN HEMISPHERE SWD03910 
JULY, ETC. UNITS ARE CM AT STP. SWD03920 

^ SWD03930 

DATA OLJAN/ . 2292 . . 230*» . . 2354 , . 241 7 . • 252 1 , . 2646 , . 2783 , . 2942 , . 3042 , SUD03940 

* . 3121 . . 3204 , . 3292 , . 3404, . 3496, . 3542 , . 3575 , . 3579 . . 3567 . . 3558/ SWD03950 

SWD03960 

DATA OLAPR/ . 2375 , . 2408. . 2475 , . 2583 , . 2725 , . 2879 , . 3062 , . 3250 , . 3429, SWD03970 

* . 3608 . . 3762 . . 3925 , . 4075 , . 4200 , . 4287 , . 4333 . . 4342 . . 4329 ..4312/ SWD03980 

SWD03990 

DATA OL JUL/ . 2387 ,. 2454 ,. 2508 , . 2583 , . 2658 , . 2 746 , . 2837 . . 2950 , . 3067 . SWD04000 

* .3187,. 3275,. 3329,. 3354,. 3358,. 3337,. 3321,. 3283,. 3229,. 3175/ SWD04010 

SWD04020 

DATA OLOCT/ . 2346 , . 2358 ,. 2383 .. 2425 ,. 2479 .. 2525 , . 2567 . . 2608 . . 2646 , SWD04030 

* .2679,. 2717 ,. 2754,. 2792,. 2829.. 2867.. 2883 . ?896, .2896, .2883/ swn04040 

SWD04050 

FOUR DISTRIBUTIONS OF TOTAL OZONE, FROM A MAX OF 0.22 TO A MAX OFSWD04060 


0.46, IN CENTIMETERS AT STP. 


DATA OCM22/ . 00008 . . 00657 , . 01830, . 03353 , . 05614 , . 08685 , . 10930 , 

* . 14029, . 16624, .17797 , . 18492 , . 18867 ,. 19120, . 19384, . 19645 , 

* .19844,. 20262,. 20601,. 20907.. 21198.. 21473,. 21728,. 21992/ 

DATA 0CM30/. 00008, .00657, .01837, .03496, .06280, .10410, . 13398, 

* . 1 752 1 , . 2 1079 . . 22947 ,. 24222 , . 2492 / . . 2541 0 , . 2591 1 . . 26396 . 

* . 26763 . . 27503 . . 28061 , . 28528 . . 28937 , . 29307 . . 29646 . . 29996/ 

DATA 0CM38/ . 00008 , . 0065 /.. 01869 , . 03923 , . 07442 , . 12224 , . 15686, 

* . 20473 , . 24695 .. 27145 ,. 29138, . 30410 ,. 31297 , . 32208, . 33065 , 

* . 33675 . . 34802 , . 35563 . . 36162 , . 36673 .. 37134, . 37559 , . 37998/ 


SWD04070 

SWD04080 

SWD04090 

SWD04100 

SWD04110 

SWD04120 

SWD04130 

SWD04140 

SWD04150 

SWD04160 

SWD04170 

SWD04180 

SWD04190 
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I 


SWD04200 


J 



DATA 0CN46/ . 00008 , . 00657 , . 01889 04170 , . 07966 , . 12844 , . 16345 , 

SWD04210 



* .21238. .25742, .28619, .31284, .33246, .34793, .36555, .38344, 

SWD04220 



* . 39675 ,. 41902 , . 43025 , . 43784 , . 4A401 , . 44958 , . 45471 , . 46000/ 

SWD04230 


c 


SHD04240 


c 

REFERENCE PRESSURE LEVELS. 

SUD04250 


c 

DATA PROCM/1. ,3. ,5 . ,7. , 10. , 15. , 20. ,30. ,45. ,60. ,80. , 100. , 120. , 150. 

SWD04:>60 

.SWD04270 



*190. ,230. ,340. ,450. ,560. ,670. .780. ,890. , 1013.25/ 

SWD04280 


c 

DATA CDATE/-77. 0,15. 0,105. 0,196. 0,288. 0,380.0/ 

SWD04290 

SWD04300 


c 

DATA TOTOZ/0. 22, 0.30, 0.38, 0.46/ 

SWD04310 

SWD04320 


c 

DATA N02/23/ 

SWD04330 

SWD04340 


c 


SWD04350 

1 

c 

note *** CALENDAR DAY-GARP REFERENCE DAY-63 *** 

SWD04360 

1 

c 

S. HEM. OZONE DISTRIBUTION IS SEASONAL REFLECTION OF N. HEM.DISTRSWD04370 

1 

c 

CDAY-XDAY-63.0 

SWD04380 

SWD0439J 

1 


IF (XLAT.lt. 0.00) CDAY-CDAY+ 183.0 

SWD04400 



IF(CDAY.GT.365.) CDAY-CDAY-365.0 

SWD04410 

1 


dlaT' abs(xlat) PAGc IS 

SWD04420 




SWD04430 

SWD04440 



IF(CDAY.LT.CDaTEJ )goto no 

SWD04450 

$ 

100 

CDATEI-CDATEJ 

SWD04460 

\ 

no 

DXDATE- (cday-cdatei) / (CDATEJ-CDATEI) 

SWD0A470 

1 


DLATI-0.0 

SWD04480 

■1 

1 


DO 120 K-2,19 

SWD04490 


J-K 

SWD0A500 

i 


DLATJ-DLATI+5.0 

SWD04510 



if (dlat.lt. DLATJ )G0T0 130 

SWD04i'20 

i 

120 

DLATI-DLATJ 

SWD04i30 

130 

DXDLAT- (DLAT-DLATl) /5 . 0 

SWJJ04540 



1-J-l 

SWD04550 


c 


SWD04560 

1 

c 

TOTAL VERTICAL OZONE CONTENT (CM*NTP) FOR GIVEN LATITUDE AND DATESWD04570 

1 

c 


SWD04580 

1 


IF(CDAY.GT. 15.00 )GOTO 150 

SWD04590 


uo 

ODl-OLOCT (1) +DXDLAT* (OLOCT (J) -OLOCT (I) ) 

SWD04600 

\ 


0D2-0LJAN (1) +DXDLAT* (OLJAN ( J) -OLJAN (I) ) 

SWD04610 

1 


GOTO 200 

SWD04620 

1 

c 


SWD04630 

i 

i 

* 

» 

150 

TF(CDAY,GT. 105.0 )G0T0 17^ 

SWDC':640 


OOnOLJAN (I) +DXDLAT* (OLJAN (J) -OLJAN (I) ) 

SWD04650 

i 


OD2 “OLAPR ( I ) +DXDLAT* (OLAPR ( J ) -OLAPR ( 1 ) ) 

SWD0A660 

1 

7 


GOTO 200 

SWD04670 

1 

■/ 

c 


SWD04680 

1 

170 

1F(CDAY.GT. 196.0 )GOTO 190 

SWD04690 



ODl-OLAPR (1) +DXDLAT* (OLAPR ( J) -OLAPR (I) ) 

SWD04700 


0D2-0LJUL(1)+DXDL.\T*(0LJUL(J)-0LJUL(I)) 

SWD04710 



GOTO 200 

SWD04720 

1 

k 

c 


SWD04730 

1 

190 

IF(CDAY.GT. 288.0 )GOTO l40 

SWD04740 



0D1-0LJUL(1)+DXDLAT* (OLJUL(J)-OLJUL(l) ) 

SWD04750 



OD2-OLOCT (I) +DXDLAT* (OLOCT (J) -OLOCT (I ) ) 

SWD04760 


c 


SWD04770 


200 

TOTOCM-ODl+DXDATE* (OD2-OD1) 

SW;JOA780 


c 


SWU04790 
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c 
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VERTICAL OZONE DISTRIBUTION FOR GIVEN TOTAL OZONE CONTENT 


210 

220 


230 

2A0 

250 

260 

270 

C 

C 

C 

280 


290 


300 


310 


C 

C*** 

C 

c 

c 

c 

c 

c 

c 

c 

c 


ORiG;r^!AL PAGE IS 
OF POOR OOALITY 


TOTOZI-TOTOZ(l) 

DO 210 K-2,4 
J-K 

TOTOZJ-TOTOZ(J) 

IFCTOTOCM.LT.TOTOZJ )G0T0 220 
TOTOZI-TOTOZJ 
I-J-1 

DXOCM- (TOTOCM-TOTOZI) /O. 08 
1F(J.GT.2 )GOTO 240 
DO 230 I-l.NOZ 

OCHXX ( 1 ) -OCH22 ( I ) +DXOCM* (OCH30 ( I ) -OCM? " ( I ) ) 
GOTO 280 

IFO.GT.3 )G0T0 260 
DO 250 1-1, NOZ 

OCMXX ( 1 ) -OCM30 ( I ) +DXOCM* (0CM38 ( 1 > -OCM30 ( Z ) ) 
GOTO 280 
DO 270 1-1, NOZ 

OCKXX ( I ) -OCM38 ( I ) +DX0CM^ ' DCM46 ( I ) -0CM38 ( I ) ) 
OZONE COi'TENT (CM*NTP) ABOVE EACH LAYER EDGE 
NP-2 

PROCMI-PROCM(l) 

DO 310 N-l.NLAYOl 
PLEN-PLE(N) 

DO 290 K-NP.NOZ 
J-K 

PROCMJ-PROCM(J) 

IF(PLEN.LT.PROCMJ )GOTO 300 

PROCMI-PROCMJ 

PROCMI-PROCM(J-l) 

DXPRO- (PLEN-PROCMI) / (PROCMJ-PROCMI) 

I-J-1 

PROCMl-“PROCM(I) 

NP-I 

OZALE (N) -OCMXX ( I ) +DXPRO* (OCMXX ( J) -OCMXX ( 1 ) ) 

RETURN 

END 

SUBROUTINE SOLAP.l (NLAY,XDAY,XLAT,XSURF') 


SHD048C0 

SWD04810 

SVID04820 

SWDO4830 

SWD04840 

SWD04850 

SHD04860 

SWDO4870 

SWD04880 

SWD04890 

SWD04900 

SWD04910 

SWD04920 

SWD04930 

SWD04940 

SWD04950 

SWD04960 

SWD04970 

SWDO4980 

SWD04990 

SMD050QO 

SWD03010 

SWD0502C 

SWD05030 

SWD05040 

SWD05050 

SWD05060 

SHD05070 

SWD05080 

SWD05090 

SWD05100 

SWDCiilO 

SWD05120 

SWD05130 

SWD05140 

SWD0S150 

SWD05160 

SWD05170 

SWD0S180 

SWD05190 

SWD05200 

SWD05210 

SWD05220 

SWD05230 


ROGER DAVIES SOLAR RADIATION PARAMETERIZATION. TAKES INTO 
ACCOUNT THE VARIATION OF SURFACE ALBEDO Hl'^d ZENITH ANGLE. 

INCLUDES ABSORPTIOr BY OZONE AND WATER VAPOR. AS WELL A3 RAYLEIGHSWD052^)C 
SCATTERING. CLOUD SCATTERING IS HANDLED IN SUBROUTINE CLOUDY. SWD05250 

SWD05260 

OUTPUT VARIABLES ARC AS(NUY), THE SOLAR FLUX ABSORBED BY LAYER SWD05270 
NUY, IN LANGLEYS PER DAY*, SG, THE SOLAR RADIATION ABSORBED BY THSWD05280 
GROUND. ALSO IN LANGLEYS PER DAY; AND TOPABS. THE SOLAR ABS0RPT1OSWD05290 
ABOVE THE MODEL TOP, ALSO IN LANGLEYS PER DAY. SWD05300 

SWD05310 

LOGICAL ICE. OCEAN. SNOW SWD05320 

SWD05330 

COMMON/RADCOM/AS(15) ,RE(16) ,PL(15) ,PLE(16) ,PLK(15) ,PLKE(16) ,TL(15) SWD05340 
*,TLE(16) ,TG,TH(15) ,SHL(15) ,SHLF.(16) , SHG, CLOUD (15) .COSZ.SO.SG.CXL, SWD05350 
* OCEAN, ICE, SNOW SWDC5360 

SWD05370 

C0MM0N/CLDC0M/SWALE(16) ,.SWIl(15) ,AL(16) .TAUL(16) ,OZALE(16) , TOPABS, SWD05 380 
1 COSHAG, SCOSZ.FSCAT.RSURF, AG SWD05390 
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c 

c 

c 

c 

c 


c 

c 


ORlGJr^Ai 

OF POOR QUALITY 


REAL TCOND( 9 ), IRENE ( 9) 

DATA TCONO/0. 0,1. 0,2.0. A. 0,6. 0,6. 0,8. 0,8. 0,8.0/ 

DATA TPEHE/0. 0,0. 0,8. 0,8. 0,8. 0,8. 0,8. 0,8. 0,8.0/ 

DATA TLOWL/16./ 

DATA TMIDL/8.0/ 

DATA NLAY02/5/ 

DATA SHLTOP , TLTOP , GTOPO , DELTA/0 . 00002 ,220.0,120.1612,0. 00001 / 
DATA PI/3.141592653/ 

GTOPO IS USED FOR SCALING WATER VAPOR ABSORPTION WITH RESPECT 
TO A STANDARD TEMPERATURE (2?3..6 K) AND PRESSURE (1013.25 MB) 

THE FORMULA FOR GTOPO IS: 

GTOPO-2 . 0*0.980*1013 . 25/SQRT (273 . 16) 

NIGHT-SIDE 

SOLAR RADIATION ABSORBED THE THE GROUND. 

SG=0.0 

SOLAR RADIATION ABSORBED BY LAYER N. 

DO 100 N=1,MLAY 
AS (N) =0.0 

1 F (COSL . LT . 0 . 01 ) RETURN 

NLAY1=NLAV^ 1 
RMEAN=0.0 

COSMAG-35 . 0/SQRT (I 2?4 . 0*C0SZ*C0S2+1 . 0) 

SURFACE REFLECTIVITY. 

AG IS THE SURFACE ALBEDO FOR DIFFUSE RADIATION. 

AG-XSURF 


RSL'RF IS THE SURFACE ALBEDO TO DIRECT RADIATION. 

IT INCLUDES THE SOLAR ZENUH ANGLE DEPENDENCE OF PALTRIDGE AND 
PLATT. 

ARG=-18.0* (PI *0,5-AC0S (COSZ) ) /P I 
AKG=AMAX1 (ARC. -75.0) 

RSURF-AG+ ( 1 . 0-AG) *EXP (ARC.) 

PARTITION OF INCIDENT FLUX SUBJECT TO SCATTERING 

SC0SZ-S0*C0SZ 

FSCAT-0.647*SCOSZ 


SWD05A00 

SWD05410 

SWD05420 

SWD05430 

SWD05440 

SWDUS450 

SWDO5460 

SWD05A70 

SWD05480 

SWD05490 

SWD05500 

SWD05510 

SWD05520 

SWD05530 

SWD05540 

SWD05550 

SWD05560 

SWD05570 

SWD05580 

SWD05590 

SWD05600 

SWD05610 

SWD05620 

SWD05630 

SWD05640 

SWD05650 

SWD05660 

SWD05670 

SWD05680 

SWD05690 

SWD05700 

SWD05710 

SWD05720 

SWD05730 

SWD05740 

SWD05750 

SWD05760 

SWD05770 

SWD05780 

SWD05790 

SWD05800 

SWD05810 

SWD05820 

SWD05830 

SWD05840 

SWD05850 

SWD05860 

SWD05870 

SWD05880 

SWD05890 

SWD059CO 

SWD05910 

SWD05920 

SWD05930 

SWD05940 

SWD05950 

SWD0596& 

SWU'J‘970 
SWI)L5S80 
Sv-O. I<990 


C 
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c 

c 


120 

C 

C 

C 

C 

C 

C 

c 


c 


c 


c 


130 


140 

C 


SCALED WATER VAPOR CONTENT ABOVE EACH LAYER EDGE 
DB-PLE(1)**2 

SWALE (1 ) -DB*SHLTOP/SQRT (TLTOP) /GTOPO 
DO 120 N-l.NLAY 
M-N+1 

DA-PLE(H)**2 

SHTEH-AHAXl (SHL (N) . 1 . 0E~8) 

W- (DA-DB) *SHTEM/SQRT (TL (N) ) /GTOPO 
SW1L(N)-W 

SWALE (M) -swale (N)+W 
DB-DA 

COMPUTE CLOUD PARAMETERS: 

TOP CLOUD, TOTAL AND FRACTIONAL CLOUDINESS PARAMETERS 

WITH N-2, CLOUDS ARE EXCLUDED FROM TOP LAYER 

FRACTIONAL CLOUDINESS IS INCLUDED, BUT NOT CURRENTLY USED. 

FCLOUD-0.0 

FCLEAR-1.0 

NTOPT-NUYl 

NTOPF-NLAYl 

DO 140 N-2,NLAY 
XX-CLOUD(N) 

IF(XX.LT.0.01 )GOTO 140 
IF(XX.GT.0.99 )G0T0 130 

FC-AMAX1(XX,FCL0UD) 

FCLOUD-FC 

FCLEAR-l.O-FCLOUD 

IF (NTOPF.lt. NLAY )GOTO 1^0 

NTOPF-N 

GOTO 140 

IF (NTOPT.lt. NLAY )GOTO 140 

NTOPT-N 

FCLEAR-0.0 

CONTINUE 

IF(FCLEAR.GT.0.99 )G0T0 200 
IF (FCLOUD.lt. 0.01) FCLOUD-1.00 


Or 




pcc: 


QUALITY 


LARGE SCALE CONDENSATION (STRATIFORM) CLOUD PARAMETERS 

DO 150 N-1,NLAY 
TAUL(N)-0.0 

IF(CLOUD(N) .GT.0.99) TAUL(N)-TC0NDC<) 

LOW LEVEL CONVECTION CLOUD PARAMETERS: 

IF(NAB.EQ.l) CLOUD IS IN UYER 7 
IF(NAB.EQ.2) CLOUD IS IN LAYER 8 

NAB-CLOUD (NLAY+ 1 ) +DELTA 
1F(NAB.LT.1 )GOTO 160 
N-NAB+6 
TAUL(N)-TLOWL 


SWD06000 

SWD06010 

SWD06020 

SWD06030 

SWD06040 

SWD060S0 

SWD06060 

SWD06070 

SWD06080 

SWD06090 

SWD06100 

SWD06110 

SWD06120 

SWD06130 

SWD06140 

SWD06150 

SWD06160 

SWD06170 

SWD06180 

SWD06190 

SWD06200 

SWD06210 

SWD06220 

SWD06230 

SWD06240 

SWD06250 

SWD06260 

SWD06270 

SWD06280 

SWD06290 

SWD06300 

SWD06310 

SWD06320 

SWD06330 

SWD06340 

SWD06350 

SWD06360 

SWD06370 

SWD06380 

SWD06390 

SWD06400 

SWD06410 

SWD06420 

SWD06430 

SWD06440 

SWD06450 

SWD06460 

SWD06470 

SWD06480 

SWD06490 

SWD06500 

SWD06510 

SWD06520 

SWD06530 

SWD06540 

SWD06550 

SWD06560 

SWD06570 

SWD06S80 

SWD06590 


c 

c 

c 

c 

160 


C 

C PENETRATING CONVECTION CLOUD PARAMETERS 

C 

170 NAC-CL0UD(NUY+3)+DELTA 
1F(NAC.LT.1 )G0T0 190 
N-NAC+3 
DO 180 1-1,4 
C 

C OPTICAL THICKNESS CALCULATION. 

C 

TAUL(N)-TPENE(N) 

180 N-N+1 

190 IF(FCLEAR.LT.0.01 )GOTO 250 
C 

C RAYLEIGH SCATTERING AND CLEAR ATMOSPHERE REFLECTIVITY 

C 

200 RBRAY-. 433/(1. 0+6. 43*C0SZ) 

RBBRAY-0.093 

RCLEAR-RBRAY+ ( 1 . 0-RBRAY) * ( 1 . 0-RBBRAY) *RSURF/ ( 1 . 0-RBBRAY*AG) 
C 

C ABSORPTION BY WATER VAPOR IN CLEAR ATMOSPHERE 
C 

W-SWALE(l)*COSMAG 
DB-AWATER (W) 

T0PABS-DB*FCLEAR*SC0S2 
DO 210 N-l.NUY 
W- SWALE (N+ 1) *COSMAG 
DA-AWATER (W) 

AL(N)-DA-DB 
210 DB-DA 

TRANS-l.O-DB 

RF-TRANS*RSURF 

AL(NLAYl) - (TRANS-0.647) * (1 . 0-RSURF) 

IF(RF.LT. 0.001 )G0T0 230 
WW-W* (1.0+1. 66/COSMAG) 

C 

DO 220 N-l.NLAY 
M-NLAYl-N 
•WW-1.66*SWALE(M) 

DA-AWATER (W) 

AL (M) -AL (M) + (DA-DB) *RF 
220 DB-DA 
C 

230 ACLEAR-FCLEAR*SCOSZ 
C 

DO 240 N-l.NLAY 
240 AS(N)-ACLEAR*AL(N) 

C 

C SOLAR FLUX ABSORBED BY THE GROUND. 

C 

SG-FCLEAR* (FSCAT* ( 1 . 0-RCLEAR) +SCObZ*AL (NLAYl) ) 

C 


MID LEVEL CONVECTION CLOUD PARAMETERS 
IF(NBC.EQ.l) CLOUD IS IN UYER 
IF(NBC.EQ.2) CLOUD IS IN LAYER 

NBC-CLOUD(NLAY+2) +DELTA 
IF(NBC.LT.l )GOTO 170 
N-NBC+4 
TAUL(N)-TM1DL 




ORiC..4.^t ^ 

OF poo^ quality 


SWD06600 
SWD06610 
SWD06620 
SWD06630 
SWD06640 
SWD06650 
SWD06660 
SWD06670 
SWD06680 
SWD06690 
SWD06700 
SWD06710 
SWD06720 
SWD06730 
SWD06740 
SWD06750 
SWD06760 
SWD06770 
SWD06780 
SWDO6790 
SWD06800 
SWD06810 
SWD06820 
SWD06830 
SWD06840 
SWD06850 
SWD06860 
SWD06870 
SWD06880 
SWD06890 
SWD06900 
SWD06910 
SWD06920 
SWD06930 
SWD06940 
SWD06950 
SWD06960 
SWD06970 
SWDO6980 
SWD06990 
SWD07000 
SWD07010 
SWD07020 
SWD07030 
SWD0704C 
SWD07050 
SWD07060 
SWDO 7 O 7 O 
SWD07080 
SWDO 7 O 9 O 
SWDO 7 IOO 
SWD07110 
SWDO 7 I 2 O 
SWD07130 
SWD07140 
SWD07150 
SWn07160 
SWOO 7 I 70 
SW007180 

swpoyno 
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c 

c 

c 

250 


260 

C 

C 

C 

C 


270 

C 

C 

C 

280 

C 

C 

C 

C 




C 

c 

c 

310 

C 


RME^N-FCLEAR*RCLEAR 
IF(FCLEAR.GT.0.99 )GOTO 290 







'c 


* fS 
JTY 


ABSORPTION BY WATER VAPOR IN CLOUDY ATHOSPHERC 

NTOP-MINO (HTOPT , NTOPF) 

CALL CLOUDY (RCLOUD.NLAY.NTOP) 

DO 260 N-l.NLAY 

AS(N)-AS(N)+AL(N) 


SOLAR FLUX ABSORBED BY THE (XOUND. 

SG-SG+FCLOUD*FSCAT* (1 . 0-RCLOUD) ♦AL (NUYl) 

RMEAN-RNEAN^FCLOUD*RCLOUD 
IF(FCL0UD.GT.0.99 )G0T0 300 
NTOP-NTOPT 
FCLOUD-l.O-FCLOUD 
DO 270 N-l.NLAY 

1 F (CLOUD (N) . LT . 0 . 99) TAUL (N) -0 . 0 
CALL CLOUDY (RCLOUD.NLAY.NTOP) 

SOLAR radiation ABSORBED BY LAYER N. 

DO 280 N-l.NLAY 
AS(N)-AS(N)+AL(N) 

SOLAR FLUX ABSORBED BY THE GROUND. 

SG-SG+FCL0UD*FSCAT* (1 .0-RCLOUD)+AL(NLAYl) 

RMEAN-RMEAN+FCLOUD*RCLOUD 
GOTO 300 


EFFECTIVE CLEAR SKY RAYLEIGH ALBEDO FOR OZONE ABSORPTION 


RBBROZ-0.144 

RBROZ-0 . 2 186/ ( 1 . CKO . 816*COSZ) 

RMEAN -RBROZ+ (1 . O-RBROZ) * (1 . O-RBBROZ) *RSURF/ (1 . 0-RBBROZ*AG) 
ROUTINE PROVIDES OZONE CN ABOVE EACH UYER EDGE 
NLAYOl-NUYOZ+1 
CALL 0Z0NE2 (NLAYOl , XDAY . XLAT) 

COMBINED UV-^VIS OZONE ABSORPTION OF INCOMING SOLAR RADIATION 

W-OZALE(l)*COSMAG 

DB-AOZONE(U) 

TOPABS-TOPABS+DB*SCOSZ 
DO 310 N-l.NLAYOZ 
W-OZALE(N+l)*COSMAG 
DA-AOZONE(W) 

DIRECT SOLAR RADIATION ABSORBED BY OZONE IN LAYER N. 

AS (N) -AS (N) + (DA-DB) *SCOSZ 
DB-DA 


SWD07200 

SWD07210 

SHD07220 

SWD07230 

SWD07240 

SWD07250 

SWD07260 

SWD07270 

SWD07280 

SWD07290 

SWD07300 

SWD07310 

SWD07320 

SWD07330 

SWD07340 

SWD07350 

SWD07360 

SWD07370 

SWD07380 

SWD07390 

SWD07400 

SWD07410 

SWD07420 

SWD07430 

SWD07440 

SWD07450 

SHD07460 

SWD07470 

SWD07480 

SWDO749O 

SWD07500 

SWD07510 

SWD07520 

SWDO753O 

SWD07S40 

SWD07550 

SWD07560 

SWDO 757 O 

SWD07580 

SWDO759O 

SWD07S00 

SWD07610 

SWD07620 

SWD07630 

SWD07640 

SWD07650 

SWD07660 

SWD07670 

SWD07680 

SWD07690 

SWD07700 

SWD07710 

SWD07720 

SHD07730 

SWD0774C 

SKD07750 

SWD07760 

SWD07770 

SHD07780 

SHDO779O 
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r 


c 

c 

c 


WW-W<»^1.90*W/COSMAG 

COMBINED UV4-VI8 OZONE ABSORPTION OF REFLECTED SOLAR RADIATION 
RF-SCOSZ*RMEAN 


DO 320 N-1,NLAY0Z 
M-NLAYOl-N 
W-WW-1.90*OZALE(M) 
DA-AOZONE(H) 


OnfGU'iaL P.’ ts 
OF roc n DUALITY 


C 

C 

C 

370 


C 

C*A* 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


REFLECTED SOLAR RADIATION ABSORBED BY OZONE IN UYER N. 

AS (M) -AS (M) ♦ (DA-DB) *RF 
DB-DA 

DA-AOZONE(WW) 

TOPABS-TOKABS+RF* (DA-DB) 

RETURN 

END 

SUBROUTINE SUNDAY (JD^Y.SOMULT.SO.DEC.SIND.COSD.EQT) 

THIS SUBROUTINE IS USED TO INITIALIZE THE ASTRONOMICAL VARIABLES 
THAT DEPEND ONLY ON THE TIME OF YEAR. 

JDAY IS THE JULIAN DAY, -1 FOR JAN 1ST. 

SOMULT MULTIPLIES THE SOLAR CONST.. TYPICALLY - UNITY. 

SO IS THE SOLAR INSOLATION AT THE TOP OF THE ATMOSPHERE, 

IN LANGLEYS PER DAY. 

DEC IS THE DECLINATION ANGLE OF THE SUN, IN RADIANS. 

SIND AND COSD ARE THE SINE AND COSINE OF DEC. 

EQT IS THE 'EQUATION OF TIME', I.E. THE OFFSET TO LOCAL TIME, IN 
HOURS. 

DATA PI/3.U1592654/ 

DTR-PI/180. 

RELATIVE POSITION OF EARTH IN ORBIT. 

PHASE-2 . *P I* (JDAY- 1) /365 . 

PHASE2-PHASE*2. 

PHaSE3-PHASE*3. 

CP-COS (PHASE) 

CP2-C0S(PHASE2) 

SP-SIN (PHASE) 

SP2-SIN(PHASE2) 

EFFECT OF EARTH-SUN DISTANCE. (PALTRIDGE AND PLATT P.57) 
ESFCTR-1.000il+0.03A221*CP+0.00128*SP-^0.000719*CP2+0.000077*SP2 
INCIDENT SOLAR RADIATION IN LY/DAY (WILLSON, 1982) 

SO-2823 . 9*S0MULT*ESFCTR 

SOLAR DECLINATION ANGLE. (PALTRIDGE AND PLATT P.63) 

DEC-O. 006918-0. 399912*CP^0.070257*SP-0.006758*CP2+0.000907*SP2 
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SWD07800 

SWDO78I0 

SWD07820 

SWD07830 

SWD07840 

SWDO78S0 

3WD07860 

SWD07870 

SWD07880 

SWD07890 

SWD07900 

SWD079I0 

SUD07920 

SWDO7930 

SWD07940 

SWD07950 

SWD07960 

SWD07970 

SWD07980 

SWD07990 

SWD08000 

SWD08C10 

SWD08020 

SWD08030 

SWD08040 

SWD08050 

SWD08060 

SWD08070 

SWD08080 

SWD08090 

SWD08100 

SWD08110 

SWD08120 

SWD08130 

SWD081A0 

SWD08150 

SWD08160 

SUD08170 

SWD08180 

SHD08190 

SWD08200 

SWD08210 

SWD08220 

SWD08230 

SWD08240 

SWD08250 

SWD08260 

SWD08270 

SWD08280 

SWD08290 

SWD08300 

SWD08310 

SWD08320 

SHD08330 

SWD083A0 

SWD08350 

SWD08360 

SWDO8370 

SWD083S0 

SWD0339C 





o o o 


1 -0.002697*COS(PHASE3)+0.001A8*SIN(PHASE3) SWD08400 

SIND-SIN(DEC) SWD08410 

COSD-COS(DEC) SWD08<»20 

SWD08430 

THE EQUATION OF TIME (HOURS) SWD08440 

SWD08450 

EQT-(0. 000075+0. 001863*CP-0.032077*SP-0.014615*CP2'0.040849*SP2) SWD08460 

1 /(DTR*15.) SWD08470 

C SWD08480 

RETURN SWO08490 

end SWD08500 
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